Charge separation effects in magnetized electron-ion plasma expansion into a vacuum 
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Charge separation effects in the expansion of magnetized relativistic electron-ion plasmas into a 
vacuum are examined using 2-1/2-dimensional particle-in-cell plasma simulations. The electrostatic 
field at the plasma surface decelerates electrons and accelerates ions. A fraction of the surface 
electrons are trapped and accelerated by the pondermotive force of the propagating electromagnetic 
pulse, a mechanism we call the DRPA (diamagnetic relativistic pulse accelerator). This charge 
separation is enhanced as the initial plasma temperature is decreased. The overall energy gain of 
the plasma particles through the expansion strongly depends on the initial plasma temperature. 
Moreover, the electrons become relatively less energized and the ions more energized as the plasma 
temperature decreases. 

PACS numbers: 



Several kinds of plasma expansions have been investi- 
gated over the past few decades. A model of free plasma 
expansion into a vacuum is interesting for studies related 
to experiments of ion jetsplQ- Structures and instabili- 
ties of plasmas expanding into a vacuum which contains a 
uniform ambient magnetic field have also been examined 
by theoretical and numerical methods. Those studies 
are relevant to laserj)lasma experiments and space 
plasma phenomenally, Also plasma injection into a 
magnetic field has been theoretically investigated^. 

Recently, Liang et al. 0, Il0| have studied a new type 
of plasma expansion, in which a magnetized relativis- 
tic plasma freely expands into a vacuum with no ex- 
ternal magnetic field. This model of expansions poten- 
tially can be applied to a variety of astrophysical phe- 
nomena, such as gamma-ray bursts[Ioj or astrophysical 
jets. When relativistic plasmas which contain a strong 
transverse magnetic field are suddenly released, an EM 
pulse forms and begins to propagate into the vacuum. Si- 
multaneously, a portion of the surface plasma is trapped 
and accelerated by the EM pulse via the pondermotive 
force. Through the expansion, the field energy of the 
EM pulse is converted into the directed kinetic energy 
of the surface particles. This new mechanism of plasma 
acceleration is called the diamagnetic relativistic pulse 
accelerator (DRPA) @. 

Electrons and ions behave differently in the expansion 
process. Because of their much smaller mass, the surface 
electrons soon outrun the ions and charge separation oc- 
curs. In this Letter, using the results of plasma full par- 
ticle simulations, we discuss charge separation effects on 
plasma energization in the relativistic magnetized plasma 
expansion into a vacuum. 

In our particle simulation code [ill], we use a two-and-a- 
half-dimensional explicit simulation scheme based on the 
particle-in-cell method for time advancing of plasma par- 
ticles. The grid separations are uniform, Ax = Az = A e , 
where A e is the electron inertial length defined by c/u pe 
(c is the speed of light and the electron plasma fre- 
quency is ojp e = ^/e 2 no/eom e ; e is elementary charge, 



no is the initial electron density, and eo is the dielec- 
tric constant of vacuum). The ratio of the electron Lar- 
mor frequency to the electron plasma frequency D, e /u> pe 
(where VL e = eBa/m e ] Bq is the background magnetic 
field) is set equal to 10 in these simulations. The ion- 
electron mass ratio is rrii/m e = 100 and the charge ratio 
is qi/e = 1. The simulation domain on the x — z plane is 
-L x /2 <x< L x /2 and -L z /2 < z < L z /2. To prevent 
large violations of Gauss's law due to numerical noise, 
Marder's method for the electric field correction[l2j is 
adopted at every time step in the code. 

We use a doubly periodic system in x and z direc- 
tions, and the system length is L x = 24QAx = 2400c/f2 e 
and L z = 120Az = 1200c/il e . Initially, the electron 
(ion) distributions are assumed to be a relativistic (non- 
relativistic) Maxwellian respectively with spatially uni- 
form temperature. The spatial distribution of the initial 
plasma has a slab form with the lengths 6Aa; x 120Az, 
and the plasma slab is located in the center of the system. 
The background magnetic field B = (0,B ,0) initially 
exists only inside the plasma. Since the system length in 
the x direction is enough large, no plasma particles and 
waves can cross the x boundaries during this calculations 
(until m e = 1000). 

We describe here the results of two simulations; one 
with initial temperatures of /csT e = ksTi — 5MeV and 
one with initial temperatures of fcsTe = ksTi = lOOfceV^. 
The qualitative evolution of the 5MeV temperature case 
is illustrated in Figure 5 of Ref. 9] . A pulse of fast elec- 
trons moves rapidly away from the ions in the ir-direction, 
carrying an imbedded electromagnetic pulse (B y > 0, 
E z < 0) with them which permits them to E x B drift. 
Behind this pulse the B y and E z are relatively weak, and 
the expanding ions lag well behind the electron pulse. 
This charge separation leads to an E x in the region be- 
tween the electron pulse and the fastest ions. 

Spatial profiles of the EM field and the plasma density 
at tVl e — 600 are displayed in Figure ^ for the lQQkeV 
temperature case. In Fig.^a), the EM pulse is propagat- 
ing in the positive x direction. The electrostatic field E x 
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FIG. 1: Simulation results for ksT — WOkeV plasma at 
tQ, e = 600. (a) Electromagnetic field (E z and B y ), (b) current 
density (J z ), (c) electrostatic field (E x ), and (d) electron and 
ion density (n e and m) as functions of x. All of the quantities 
are obtained at z — 0. Phase plots p x /m e c— x for (e) electrons 
and (f) ions. These profiles should be compared with those 
of the k B T = 5MeV case in Ref.[8l. 



shown in Fig. [Jc) is generated by the charge separation 
between electrons and ions. As shown in Fig.^d), elec- 
trons are distributed spatially more extensively than ions 
because of the larger mobility of electrons. In the case 
of electron-positron expansion [j| , there is no difference 
in the density distributions of electrons and positrons. 
Consequently, both electrons and positrons are equally 
energized by the DRPA. But in the present electron-ion 
case, only electrons can follow and get energized by the 
EM pulse, while the ions are accelerated in a secondary 
manner by the charge-separation electric field. 

Phase plots p x /m e c — x of (e) electrons and (f) ions 
at tVL e = 600 in the \00keV temperature case are also 
shown in Fig. ^ In Fig-d e ), electrons are decelerated 
in the region 300 < x/(c/Cl e ) < 600, where the electric 
field E x exists. The surface electrons in the vicinity of 
x ~ 600c/il e are strongly energized by the EM pulse. 
Since some of the expanding electrons are reflected at the 
potential wall of the electrostatic field at x ~ 300c/f2 e , 
a loop structure in phase space can be seen. On the 
contrary, the front ions in this region are accelerated by 
E x and the averaged kinetic energy of the front ions tends 
to increase with time. The front ions are not affected by 
the DRPA since the ions cannot follow the EM pulse. 
Hence, we obtain an interesting result regarding particle 
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FIG. 2: System-integrated energy in the magnetic field, elec- 
tric field, electrons, and ions for (a) ksT = 5MeV and (b) 
fcflT = lOOfceV" plasmas. Each energy is normalized by the 
initial magnetic field energy in each figure. 



energization in the electron-ion plasma expansion: the 
surface electrons expanding in phase with the EM pulse 
are accelerated by the DRPA, but a portion of the front 
electrons are also decelerated by the electrostatic field 
E x caused by the charge separation, while the front ions 
are always accelerated by E x . Here, we showed only the 
lOOfceV temperature case in Fig. See Ref.0 for the 
profiles of the 5MeV case. 

Temporal evolutions of each energy component 
through the expansion are shown in Figure |3 for two 
cases with (a)k B T e = k B Ti = 5MeV and {b)k B T e = 
k B Ti = lOOfceT^. Note that in case (a) the initial plasma 
f3 (=plasma pressure/magnetic pressure)^ 0.216 while 
in case (b) j3 = 0.02, so the plasma is strongly magnetic- 
dominated. During < til e < 80, roughly correspond- 
ing to the light crossing time of the initial plasma slab, 
the system evolves from the initial magnetostatic config- 
uration into two counter-propagating EM pulses loaded 
with plasma. In this early phase, the magnetic energy 
decreases and the electric energy increases until the two 
are almost equal. After til e ~ 80, the field energy de- 
cays and is converted into the directed kinetic energy 
of plasma particles. In Figs. 13a) and (b), we observe 
the monotonic increase of the total kinetic energy of the 
ions at all times. The total electron energy, however, de- 
creases from tfl e = 80 to 400 in the k B T — 5MeV case, 
and hardly changes after tVL e = 400 in Fig. Ufa). Even 
for the lOOkeV case of Fig. Efb), though the total elec- 
tron energy increases monotonically, we see that the rate 
of increase of electron energy is much smaller than that 
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FIG. 3: Temporal evolution of peak value of electrostatic field 
E x . Comparison between different temperature plasmas. 



of the ion energy. In addition, the decrease of the total 
electron kinetic energy in Fig. |2Ia) is mainly caused by 
the deceleration due to the electrostatic field E x (we de- 
scribed the details in Figure^]). The same thing happens 
also in Fig. |5Jb). But because the initial thermal energy 
of electrons is less than the directed kinetic energy given 
by the DRPA, the overall kinetic energy of electrons in- 
creases monotonically for the lOOkeV plasma. 

In Figure |21 the peak value of E x for different initial 
temperatures are shown as a function of time. The elec- 
trostatic field E x peaks in the vicinity of the ion front. 
At the early phase of the plasma expansion (if2 e ~ 180 
in Fig. |3J) , the field E x is formed by the charge separa- 
tion. After that period, the density difference between 
electrons and ions is gradually reduced by ion accelera- 
tion and electron deceleration due to E x . As a result of 
this, the field E x becomes weaker. A plasma tends to be 
spatially broader as its temperature becomes higher. Ac- 
cordingly, the higher the plasma temperature, the smaller 
the density difference between electron and ion becomes. 
That is why the peak value of E x in the case of the 
\QQkeV temperature plasma tends to be larger than that 
of the 5MeV temperature plasma in Fig. |21 

Let us evaluate quantitatively how the field E x influ- 
ences particle energization. Fig. ^ displays the momen- 
tum distribution of the front electrons at three differ- 
ent times for initial plasma temperatures of (a) ksT — 
hMeV and (b) k B T = lOOkeV. The front electrons are 
accelerated by the DRPA and simultaneously decelerated 
by the electrostatic field as the expansion proceeds. The 
energy gain of these electrons seems to be positive and 
the kinetic energy continues to increase as time elapses, 
as shown in Fig. 

This means that the acceleration process dominates 
the deceleration process overall for the front electrons. 
Figs. 0Ja) and^Ib) show that the electron energy gain 
depends on the initial plasma temperature. The electric 
field induced by charge separation tends to be smaller 
when the plasma temperature is higher. Thus, the front 
electrons are more strongly decelerated in the low tem- 
perature case, and the energy gain of the electrons be- 
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FIG. 4: Momentum distributions of electrons in the ex- 
pansion front p x /m e c — x for (a) k B T — 5MeV and (b) 
k B T = 100/fceV plasmas at tSl e = 0, 400, and 1000. 
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FIG. 5: Momentum distributions of ions in the expansion 
front p x /m e c-x for (a) k B T = 5MeV and (b) k B T = WOkeV 
plasmas at ifi e = 0, 400, and 1000. 



comes large when the temperature is high. 

Figure displays the momentum distribution of the 
ions in the expansion front for (a) ksT = bMeV and (b) 
ksT = lOOfceT^ plasmas. Contrary to the case of elec- 
trons, the peak value of the distribution becomes smaller 
as the initial temperature increases [The peak values are 
p x ~ 140m e c in Fig.[5fa) and p x ~ 170m e c in Fig. Efb), 
respectively] . The front ions are accelerated only by the 
electrostatic field of the charge separation. The energy 
gain of the front ions tends to be large when the plasma 
temperature is low since the electric field becomes small 



4 



as the temperature increases. 

In summary, we have examined charge separation ef- 
fects in expanding electron-ion plasmas in the context 
of the diamagnetic relativistic pulse accelerator (DRPA) . 
When the magnetized electron-ion plasma expands into a 
vacuum with no external magnetic field, only the surface 
electrons are efficiently accelerated by the DRPA. Be- 
cause of the difference of the mobility between electrons 
and ions, charge separation occurs and generates a strong 
electrostatic field. This electric field accelerates the front 
ions and decelerates some of the front electrons. The 
initial plasma temperature affects the charge separation 
effect. As the temperature increases, the electric field 
caused by the charge separation becomes small. There- 
fore, electrons in high temperature plasmas are energized 
by the DRPA more efficiently compared to electrons in 



low temperature plasmas. On the contrary, ions are en- 
ergized by charge separation more efficiently in the low 
temperature case. 
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